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We have computed Bk, using two different methods with staggered fermions on a 16'^ x 40 lattice at /3 = 5.7 
with two dynamical flavors of a mass 0.01. Using an improved wall source method, we have studied a series of 
non-degenerate quark antiquark pairs and observed no effect on Bk, although effects were seen on the individual 
terms making up Bk- 



1. INTRODUCTION 

The standard model which is believed to de- 
scribe our world of hadrons and leptons contains 
a number of fundamental parameters. The knowl- 
edge of hadronic weak matrix elements is crucial 
to determine these parameters from the experi- 
ments. Lattice gauge theory has reached a point 
that it is capable of calculating hadronic weak 
matrix elements. Especially, the knowledge of 
Bfc which describes the neutral K meson mix- 
ing is crucial to narrow the domain of Vtd and 
the top quark mass, which are the fundamental 
parameters of the standard model. 

We have computed Bk on a IG'^ x 40 lattice 
at ^ = 5.7 (a"i ~ 2.0 GeV) with two dynami- 
cal flavors of a mass 0.01. The results were ob- 
tained over 155 gauge configurations. Our work 
extends earlier calculations of Bk and includes 
experiments with alternative lattice formulations 
and improved wall source methods. There are 
two methods to transcribe the continuum weak 
matrix elements to the lattice with staggered 
fermions I^J^ : the one spin trace formalism and 
the two spin trace formalism. The two spin 
trace formalism (2TR) has been used predomi- 
nantly for weak matrix element calculations 
Recently, the one spin trace formalism (ITR) has 
been developed to a level which allows it to be 
used for weak matrix element calculation on the 



lattice [|l|. We have tried both formalisms to 
calculate Bk- The results are compared in this 
paper. We have studied an improved wall source 
named cubic wall source in order to project out a 
specific hadronic state exclusively The results 
are compared with those of the conventional even- 
odd wall source. We have looked into the effect 
of non-degenerate quark antiquark pairs on Bk 
and the individual components making up Bk in 
detail. From the standpoint of chiral perturba- 
tion theory, the effects of non-degenerate valence 
quark are related to the rj' hairpin diagram in 
(partially) quenched QCdQ. We interpret our 
numerical results in terms of chiral perturbation 
theory. 

Preliminary results have already appeared in 
Ref. §]. 

2. K - K Mixing 

In the continuum, Bk is defined as 
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The operator transcription of the continuum 
AS* = 2 operator to the lattice in both one spin 
trace formalism and two spin trace formalism is 
explained in Ref. jl]. The numerical results of 
Bk calculated in both formalisms are compared 
in Figure |^. The details of renormalization proce- 
dure is given in Ref. |^] . From Figure |l|, note that 
numerical results of the tadpole-improved renor- 
malized Bk in both formalisms agree with each 
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Figure 1. Comparison of Bk in one spin trace 
form (filled circle) with Bk in two spin trace form 
(empty square) . Top: unrenormalized Bj^ . Bot- 
tom: tadpole- improved renormalized Bk at tt/u 
scale. Even-odd wall source is used. 
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Figure 2. Unrenormalized Bk of wrong (scalar- 
like) flavor structure {{V + A) ^ S')^™ and of 
a valence quark mass 0.02 versus the Euclidean 
time. Calculated in two spin trace form. Top: 
even-odd wall source is used. Bottom: cubic 
wall source is used. 



other better than those of unrenormalized Bk- 

Here, we address two important questions on 
the validity of our approach to Bk- The 
higher loop radiative correction of four fermion 
operators cause the mixing of operators with dif- 
ferent spin and flavor structures. The first ques- 
tion is how much the mixing of operators with 
different flavor structure contributes to our weak 
matrix measurement at finite lattice spacing non- 
perturbatively. The second question is how exclu- 
sively we can select the pseudo-Goldstone mode 
using our improved wall source named cubic wall 
source- We have chosen {{V + A) ® S)"^"^^' in order 
to address the above two questions, where we fol- 
low the notation in Rcf. [p|. The matrix element 



of this operator with K mesons is supposed to 
vanish in the continuum limit (a 0) of lattice 
QCD, due to vanishing fiavor trace. We present 
the numerical results of the wrong fiavor channel 
in Figure ||, which tells us that the wrong flavor 
channel is highly suppressed (within 1% of Bk) 
when cubic source is used. This implies that un- 
wanted mixing of {{V + A) ^ S)^"^^ should be at 
most 1% of Bk since it is suppressed by as/(47r) 
as well as by vanishing flavor trace. The small de- 
viation of the wrong flavor channel from zero in 
Figure |^ tells us how efficiently the improved wall 
sources suppress contaminations from unwanted 
hadronic states. In Figure 0, we compare Bk 
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Figure 3. Unrenormalized Bk vs. average quark 
mass {mqo): The filled circle (empty square) rep- 
resents the data of cubic wall source (even-odd 
wall source). Calculated in the two spin trace 
formalism. 



of cubic wall source with Bk of even-odd wall 
source. 

Let us jump into the next issue: the effects 
of non-degenerate quark antiquark pairs on Bk 
and its individual components, and their rela- 
tionship with various chiral logarithms. In our 
numerical simulation, we have used three degen- 
erate quark antiquark pairs: {0.01, 0.02, 0.03}, 
and four non-degenerate pairs: {(0.004, 0.01), 
(0.04,0.02), (0.01, 0.03), (0.01,0.03)} to produce 
K meson state. In order to discuss individual 
components of Bk in an organized way, we need 
to consider a theory with four valence flavors 0: 
S and S' both with mass TOj as well as D and D' 
both with mass m^. Let be the S^^D pion 
and K'^ be the corresponding state with primed 
quarks (S'^^D'). We define individual compo- 
nents of Bk as 
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Figure 4. Unrenormalized Bk, By, Ba vs. av- 
erage quark mass: Calculated using cubic wall 
source in the two spin trace formalism. Ba- 
empty triangle and filled square. Bk- star and 
filled circle. By- x and filled diamond. Filled 
square, circle, and diamond represent degenerate 
quark antiquark pairs. All other symbols corre- 
spond to non-degenerate pairs. 



where a, b represent color indices. Let us sum- 
marize chiral perturbation results for the above 
individual components, the details of which are 
explained in Ref. ||^. In full QCD chiral per- 
turbation, the one- loop corrections to Byi, By 2, 
Bai, Sind Ba2 include a logarithmically divergent 
term in the chiral limit — * 0, which is called 
an enhanced chiral logarithms, and which is ab- 
sent in Bk- The enhanced chiral logarithms are 
not a function of quark mass difference between s 
and d quarks. We present Bk, By, and Ba data 
versus quark mass in Figure which illustrates 
the existence of a divergence in By and Ba and 
no divergence in Bk in the chiral limit. 

There are two possible differences between par- 
tially quenched or quenched and full QCD. The 
first difference is that the meson spectrum eigen- 
states are not same. The second difference comes 
from rj' loops, which can not contribute in full 
QCD because it is too heavy for chiral dynam- 
ics. The hairpin diagram contribution to Bk and 
its individual component, which is called ^^(par- 
tially) quenched chiral log" , vanishes in the limit 
of = md- In other words, (partially) quenched 
chiral logarithms arc functions of quark mass dif- 
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Figure 5. ABai vs. average quark mass: Calcu- 
lated using cubic wall source in the two spin trace 
formalism. 
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Figure 6. ABk versus average quark mass; Cal- 
culated using cubic wall source in the two spin 
trace formalism. 



ference. Since chiral perturbation theory predicts 
that Bai has a constant term about 1/3 smaller 
and a enhanced chiral logarithmic term around 
3 times larger than Ba2, we have chosen Bai 
as a useful measurement adequate to observe, if 
present, both enhanced and quenched chiral log- 
arithms. In order to observe the effect of non- 
degenerate quark antiquark pairs on Bai, hi'st we 
fit data of degenerate pairs to obtain B'^^irriK) 
and second we subtract those degenerate contri- 
butions from non-degenerate data as follows: 



ABAiimK,e) 



BAiirriK, e) - B'^^{mK 



where BAiirriKje) is data of non-degenerate 
quark antiquark pairs and e = (m^ — md)/ {rus + 
nid)- We present ABai data in Figure ^, which 
indicates that there is an additional divergence as 
a function of e. This additional divergence might 
come from (partially) quenched chiral logarithms 
or from a finite volume effect. This needs more 
numerical evidence and further theoretical under- 
standing. 

In order to detect the effect of non-degenerate 
quark antiquarks on Bk, we follow the procedure 
similar to Bai case. First, we fit degenerate date 
to B'^^irriK)- Second, we introduce ABk as fol- 
lows: 



ABK{mK) 



perturbation predicts that leading effect of non- 
degenerate quark antiquark pairs on Bk is of or- 
der e^. We present the ABk data in Figure ^, 
which illustrates that the non-degenerate quark 
mass effect is much smaller than our statistical 
error. 

Our best results are unrenormalized BximK) 
= 0.658(77) and tadpole-improved renormalized 
(N.D.R.) BKimK,tJ.= f ) = 0.659(63). 
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where the function is normalized by since chiral 



